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1 . INTRODUCTION 


The time-dependent compressible Euler equations are a hyperbolic system 
of five conservation laws 

3 

(1.1) u + l (f (u)) = 0, uf r. 

C i=l 1 X i 

+ 2 

In spherical coordinates (r,0)€R x S , the equations assume the following 
form 

1 2 1 

(1.2) U + F(e,u) + -i- l G,(0,u) +^H(0,u)-O. 

t r r i=l 1 0 i r 

Since they are rotationally invariant, there exists K(0) such that 

K _1 F(0,Ku) = F(u) 

K _1 G 1 (0,Ku) = G 1 (u) 

- I K" 1 G (0,Ku) + K -1 (0)H(0,Ku) = H(u) 

°i 

for some V , 2^, H. Then, by setting u = K(0)v we can express (1.2) as 
d.3) v t + (?(v)) r + i l e i (v) e + i H(v) = 0. 

For convenience we suppress the tildas in (1.3). A similarity solution of 

(1.3) depends only on s = — , 0 and satisfies 


(1.4) 


(F(v) - sv) s + I (G(v)) 0 + I H(v) + v = 0 
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where 


or, in nonconservative form 


G 0 " E (G i ) 0 1 ’ 


(F (v) - slv + - G (v)v. + — H(v) = 0. 
v' ' S S V 0 s 


Let L(v)F v (v) = A(v)L(v) with 


A(v) = 


AjCv) 0 


0 A^(v)j 


If we multiply by L(v) on the left 


(1.5) 


(A - s)Lv + - LG v + i LH = 0. 
s s v 0 s 


Suppose one substitutes 


( 1 . 6 ) 


v(0,s) = Vq + ev^O.s) + e v 2 + 


into (1.5), with vg a constant solution. Then Vj(0,s) will satisfy the 
linear equation 


(A(v q ) - s)L(v 0 )(v 1 ) s + j L ( V O )V V O )(V 1 ) 0 + 7 (LH) v (v 0 )v l = °* 


Suppose we had the appropriate boundary condition, and subsequently 
solved the boundary value problem for v^. Note that the system degenerates 
where s = Aj(vq), j = !,•••, 5; therefore, one suspects the solution 
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v. (0,s) will develop a singularity at that point. Later we will see that the 


singularity is typically of the form /s - X. . This shows that blows 

^ s 

up at Aj and therefore the expansion (1.6) is no longer valid when s is 

2 

too close to Aj. In our case "too close" means within e . 

In order to see nonlinear phenomena such as shocks and expansion waves, 
one needs to let s vary near A^. Thus we are led to an ansatz different 
than (1.6). 

Near an eigenvalue A^.(Vq) we rescale 


2 - 


and we let 
(1.7a) 


s = Aj(v q ) + eVAj(v o )«v 1 (0 , Aj(vq)) + e s 
v “ v Q + eVj^Q, Aj ( vq)) + e 2 v(6,s) + 


Note that if 

v^e.s) = (e , Aj (v Q ) ) + /s - Xj(v Q ) v^e) + 
then in the expansion (1.6) 


(1.7b) v = v 0 + ev^Q.A^Vq)) + e/s - A..(v Q ) v^e) + ••«. 

Our application is the shock-wedge diffraction problem in two dimensions 
(Figure 1.2), We will substitute (1.7a) and solve for "v explicitly by 
matching its boundary values with the third term in (1.7b), which is obtained 
from Keller and Blank [1] where the solution to the linearized problem is 
given. The term "v will capture the position of the diffracted shock, within 
the order of approximation. It will contain vorticity generated by the curved 
shock. The 0(e) approximation, v^ , is too crude to see such effects. 
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To be specific we consider the isentropic, two-dimensional Euler 
equations in polar coordinates and self similar form (see (1.4)): 



and the equation of state is p - Ap^. 

The assumption of constant entropy is no loss of generality because the 
changes in entropy are of order lower than the order of our approximation. 
Differentiating (1.8) and simplifying 



(1.9a) 
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Multiplying on the left by the eigenvector matrix 

P o 

-p 0 , 

0 1 




(1.9b) 



Consider now a weak shock impinging on a wedge so that at time t = 0 
the shock reaches the wedge. The initial values at t = 0, comprising of the 
states ahead and behind the shock, are shown in Figure 1.1. 
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( 1 ) 

(1) 

p = po+Ap 

u (1) = AU 

( 1 ) iA 
v = po+Ap 

( 1 ) 

c =co+Ac 


( 0 ) 

P= po 
u = o 
v= o 

P= po 

C= CO 



Figure 1.1 


The problem above could be viewed as the simplest Riemann Problem in two 
dimensions that cannot be reduced to one dimension. A Riemann Problem is an 
initial value problem with piecewise constant states and straight 
discontinuity interfaces. The existence of a solution remains an open 
question; however, it is presumed the solutions are of similarity type. We 
wish to obtain the first few terms in the uniform formal asymptotic expansion 
approximating the similarity solution in Figure 1.1. The solution is 
symmetric about the real axis. Its upper half plane restriction is shown in 
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similarity coordinates in Figure 1.2. It consists of the reflected shock 
(S R ) which interacts with the expansion wave E and curves into the 
diffracted shock Sq. 

The state (2) behind the reflected shock is constant and could be 
obtained exactly, together with the shock angle, from shock polars [4], We 
only need the first two asymptotic terras which we derive in Section 2. The 
curve separating state (2) from the expansion wave is characteristic, and the 
solution is expected to have a gradient discontinuity across it, in analogy 
with the one-dimensional rarefactions. 



Figure 1.2 
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Away from the "corner" where E meets Sq, we scale the radial variable 
in accordance with (1.7a) and obtain the asymptotic behavior of Sq and E. 
The region where E meets S D requires a special ansatz where one needs to 
scale 0 appropriately as well. This is to be expected since changes in the 
tangential direction are substantial there. As a result we obtain a transonic 
equation, resembling the small disturbance equation in steady flows, for which 
we can provide boundary condition by matching with the outside. Solving this 
boundary value problem analytically, or even proving the existence of a 
solution, remains an open question. One encounters difficulties analogous to 
the boundary value problem for the small disturbance equation. 

Hunter and Keller have applied weakly nonlinear geometrical optics, [2], 
to the shock-wedge diffraction problem [3], Their ansatz makes use of 
additional fast varying variables in analogy with the geometric optics 
construction for highly oscillatory solutions. This method, however, did not 
apply to the interaction at the "corner" (see Figure 1.2). In contrast with 
[3] we consider the equivalent boundary value problem for the pseudosteady 
solution in similarity coordinates and systematically match asymptotic 
expansions in which the original variables are rescaled to reflect nonlinear 
behavior. This way we can obtain the uniform asymptotic behavior over the 
whole space. 

We mention that in one space dimension DiPerna and Majda [5] have 
justified the weakly nonlinear geometric optics construction. 
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2. SHOCK REFLECTION 

Here we will obtain the first two terms in the asymptotic approximation 

to the shock reflection at the wall (see Figure 1.2). The small parameter, 

throughout this paper, is = e . From normal shock relations, one easily 

p 0 

computes that for state (1): 


( 1 ) 

£ 

p 0 

c 0 


= 1 + 


+ £ 2 + 0(e 3 ) 


( 2 . 1 ) 


v“> = 0 


( 1 ) 


= 1 + e + (y - ~ ^ ~ 3) e 2 + 0( £ 3 ) 


1 + xil e + 0 ( e 2 ) 

0 


where Sj- is the speed of the incident shock. 

We will now show how to obtain the asymptotic terms for state (2). The 
first term is simple and easy to compute. It corresponds to the linear theory 
of reflection. In particular it shows that the angle of reflection equals the 
angle of incidence. The second term is rather tedious to compute. 

We rotate the plane so that the wall becomes horizontal (see Figure 2.1). 

The Rankine-Hugoniot conditions across the reflected shock in Cartesian 


coordinates are: 



H 


sin g putu - x_; + pi + cos g puv 


where Xg is the velocity of the incident shock along the ramp, g is the 
angle of reflection and the subscript at the brackets denotes a jump at the 
reflected shock. 


Figure 2.1 


From (2.1) we obtain the rotated velocities for state (1) 


u v / v — Q 

= (e + — r — e ) cos 0 + 0(e ), 


= ~(e + — -T- - 3 - e 2 )sin a_ + 0(e 3 ), 
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and 


- _L_ [1 + IdJ. e) + 0(e 2 ). 

Cq cos a L 4 J 


For state (2) we set 


( 2 ) 


= 1 + ep| 2) + e 2 p^ 2) + 0(e 3 ), 


and 


!^-l +E u< 2) + E 2 u< 2 > + 0( e 3 ), 

c o 1 2 


v (2) = 0, 


tan 0 = tan g Q + e tan + 0(e ). 


Substituting in (2.2) we obtain the 0(e) equations 


2 ( 2 ) ( 2 ) 

sin g Q (-sin a Q - cos a Q u^ + p| ) - cos g Q sin ct Q cos a Q 


(2) (2) n 

u^ - cos Qq Pj^ = 0 


( 2 ) 

sin 3 q tan cig + cos gg(l - pj ) = 0 


( 2 ) - ? 

Pi “ 2 


( 2 ) 

u, =2 cos a, 


0 


and the solution 
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~ 2 V 


After some algebra we obtain the following equations from 0(e) terms: 


( 2 ) 

^2 ( 2 ) 2 1 

- u 2 + tan a 0 sin a Q tan ^ + cos a Q 

( 2 ) 

_ ( 2 ) . U 2 

2 cos a. 


= 0, 


+ 2 sin 2 a Q + 5 = 0, 


( 2 ) 2 2 

-p 2 + tan Oq tan g^ + sin Oq - cos Qq + 2 - y = 0 


( 2 ) ( 2 ) 

from which we can obtain the solution p 2 , u£ g^, as long as 0 < <*q < 

For simplicity of notation we won't write them down explicitly. In the 
original, unrotated coordinates: 


( 2 ) 


1 4 . 9 + ( 2 ) 2 

1 + 2e + $2 e » 


(2.3) 


u 


.( 2 ) 


'0 

'0 

,( 2 ) 

“0 


2 ( 2 ) 2 

2 cos e + cos u 2 e , 


94 , ( 2 ) 2 
2 sin Oq cos Oq e - sin ciq u^ e , 


= 1 + (y - l)e + ( - ~ 2 ~~ p 2 2 ^ + V 2 Cy ~ D(y ~ 3))e 2 , 


tan g = Cq + Oq + e tan g^ + 0(e ), 


g = angle of Sp with wall, and from (2.1) and (2.3) 


t'oM 
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(i) 


(i) (i) 2 

1 + e + p 2 e + 


( 2 . 4 ) 


R 


(i) 


0 

(i) 


0 


'0 


R^e + R 2 1 ^ e 2 + 


(i) . (i) 2 , 

©1 e + ©2 e + 


i = 1,2 


and 


(i) _ 


1 i = 1 


2 i = 2 


(i) _ 


( 2 ) 


i = 1 


1 = 2 


R 


(i) 


cos e 


i = 1 


2 cos dp cos(6 - Oq) i = 2 


Y " 3 


cos 0 1=1 


R 


(i) J 


R2 2) (0) 


1 = 2 


,(i) = 


-sin 0 


i = 1 


, 0 


(i) 


-2 cos c*q sln(0 - Oq) i = 2 


_ Y - 3 


©2 2 ) (e) 


sin 0 i = 1 


1 = 2 


( 2 ) ( 2 ) ( 2 ) 

where , R^ » ®2 can computed, 
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3. INTERIOR ANSA.TZ 

We start with an "interior" approximation analogous to (1.6): 


r — 
c 0 


2 

— — = 1 + eo + e o + • • • 
p 0 

(3.1) 

— — = eR + R + ••• 

C 0 


- 60 + 0 + •••. 

C 0 


Substituting in (1.9a) and collecting 0(e) terms, we obtain the linear 
equations 

-r p + R + I (0 + R) = 0 
IT IT — 9 


(3.2) p_-rR_=Q 

r r 


“ r ©_ + = P e - o, 
r r 


which imply the following equation in p alone 


(3.3) r 2 [(l - r 2 )p_]_ + p + r p_ = 0. 

r r 00 r 


Remark : Solutions to (3.3) are just similarity solutions to the wave 

equation p = Ap. 
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We go further and collect 0(e) terms and obtain 


(3.4) 


r 2 [(l “ r 2 )p_]_ + p + r p_ = G, 
r r r 


where G has nonlinear dependence on p , p_, p, p . 

r r r 0 

In order to obtain the 0(e) uniform approximation we seek solutions to 

(3.3) with ifi- = 0 on the solid body and by Huygen's Principle (see [1] and 
3n 

the above Remark) and (2.4): 


(3.5) 


p(l,0) = P 1 = 


1 2Uq < 0 IT 


2 Oq < 0 < 2<Xq 


The solution was obtained by Keller and Blank [1] using the Buseraann 
t rans f ormat i on 


1 + A - r 2 


which takes (3.3) into the Laplace equation 


ip - hUp P 66 - 0- 


The solution they obtained is 


1 + - tg" 1 


(3.6) 


(1 - & 2X )cos Xtt 


11 -(1 + £. X )sin Xtt - 2& X sin X(0 - r) 


L 1 “1 
+ - tg 


-(1 - l A )cos Xir 


77 1(1 + Jl 2X )sin Xtt - 2£ X sin X(0 - tt) 
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with X = j ( — ), and tg -1 : 3R + [0 , tt ] - Note that sin x(ir - 2 a ) 

^ '■TT CXq' U 

= sin Xir. We are interested in the asymptotic behavior of (3.6) for £ near 

1. 

For 0 * 2<Xq 

(3.7) 7=p ' / ~ X — L — - 2Xtf - / - 7 +0(1-7). 

17 (sin At:) - (sin x(r - 0)) 

The approximation (3.7) is not valid for 0 near 2ctg. We observe, however, 
that if r = 1 + er' and 0 = 2ag + e ^ 0"» then 

(3.8) “ = 1 + I tg _1 + 0 ( e 1/ 2 ), r ' < 0. 

t: 0 

Since -r^ 0 + p = 0 we also have 

- 0 

0 = 0j(2ciq) + 0(e ^ + •••• 

The two expansions (3.7), (3.8) match in the region where 0 = 2ag + p0 and 

e 1/2 « n « 1. 

Using (3.6) in (3.4) and the Busemann variables, 


Ap = G 


Gq(0 ) 

and G = j + G (0 , £) for (0,£) near (0 Q ,1) with a fixed e Q * 2a Q 

and (1 - £) Gj(0,£) continuous at £ = 1. Let p = p ~ V 2 ^ • Then 

~ -2 -2 

Ao € H, where H, is the local Sobolev space for a neighborhood of 
K loc loc 

(0 O ,D. 
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Using local regularity results for A we obtain p € Lr and therefore 

ioc 

the asymptotic behavior of p for £ near 1,0* 2ctg: 


P ~ V 2 ( 0 ) y-~ ~ V 2 G o (0) - - ■ 

A - r 


This shows that in the expansion (3.1), we have 


2 = 

p/pg =l+ep+e p+ ••• = 1 + e 


P x + eOn/l - r)j 


(3.9) 


+ E 2 0 


A - r ; 


+ ••• for 0 * 2cxq. 


In order to justify the boundary values in (3.5), that is, match (3.1) and 

2 1 

(2.4) to 0(e), it is necessary that e « e. Therefore, we restrict 

yTT? 

the linear solution (3.6) to 


(3.10) 


e « 1 - r « 1. 


We have not yet completed the 0(e) asymptotics since we haven't matched the 
"corner" expansion (3.8). We postpone it until the last section. With 0 
not near 2 (Xq we wish to move in closer to the eigenvalue and match the 
A - r term in the expansion of p (3.7). 


4. EXTERIOR ANSATZ 

In this section we consider 0 away from 2ctg. Following the procedure 

outlined in the Introduction we rescale the radial variable close to an 

2 

eigenvalue. The eigenvalue in question is R + c = Cq( 1 + c^e + R^e) + 0(e ). 
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Consequent ly our ansatz is 


r , . . 2 — 

— = 1 + c, e + R.(0)e + e r 

c 0 1 1 


— = 1 + ep. + e 2 p(r,0) + 0(e 3 ) 
p 0 1 


(4.1) — = eR.(e) + e 2 R(r,0) + 0(e 3 ) 

C 0 1 


— = £0,(9) + e 2 0(r , 8) + e 3 0(r,0) + 0(e 4 ) 
c 0 


— — = 1 + C. £ + C£ 2 + ••• 
C 0 1 


with p^, Rp 0^ from (2.4). 

Substituting (4.1) in (1.9b) we obtain the following equations in 
nondimensional (barred) variables 


e 2 (R - r + c) + 


0 

0 


eOj + 


£0 j 


e 0 ^ + ... — £0 ^ + • 

1 + • • • 0 


“ 


** “ 

0 0 


p + R + • • • 

2 + • • » 0 


p — R + • • • 

• 

• 

+ 

l 

o 


. 0 + £0 + •••- 


r 


1 + e(c^ + p^) + • • • 

( 

* 2 - 

e p + . . . 

1 + • • • 


eRj + ••• 

£0j + ••• 

( 

_£0j + £ 2 0 + • • •_ 


ep 


- R, 


eR + • • • 


£0 + e 2 0 + 
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0 

1 +e(c^ + p^) + ••• 

— £0^ + ••» 


1 + 

ep^ + * • • 

+ 

0 

1 + ... 

£0^ + ••• 


eR^ 

+ R + « • » 


_0 

£0^ + • • • 

0 


_e0j 

+ . . . 

Consider the 

third equation first. 

One obtains 

from 

0(1) terms: 


(4.2a) — 0 = 0 0 = 02 ( 0 ), 

r 

and from 0 (e) terms 


(4.2b) -0_ - Rj p_ = 0 0_ = -Rj p_. 

r 8 r r 0 r 


The second equation yields 


(4.3) 


p_ - R_ = 0 R - R2(0) = p - P 2 
r r 


from 0(1) terms. 

Since 0, + R. =0 (see (2.4)) and 0_ = 0, the first equation has 

0 r 

2 2 

only terms of 0(e) or higher. From 0(e ) terms one obtains 


(R - r + c)(p + R) + 0, R, + 0 - 0. R, (p + R) 

7 1 *0 0 1 *0 7 


- R 0_ + R - 0? = 0. 
X 0 r 
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In view of (4.2b), (4.3) and the fact that 


R i - 0* - R ± + 
6 


0, =0, i = 1,2, 


c 


Y - 1 — (y - l)(y - 3) 
2 p + 8 


one obtains 

(4.4) ((y + l)p - 2x)p + p = 0 

X 


where 



The solution to (4.4) is given by 


(4.5) 


x-(y+1)p= Kp^ 


with K = K(0) the constant of integration, from which 


(4.6) 


~ — (y + 1) ± Ay + l)^ + 4Kx 

p =— 


The solution p is uniquely determined by boundary conditions at x ±<» 
which are obtained from matching the expansion in (4.1) with (2.4) and 
(3.1). To match (2.4), p = 0 as x +<*>. This immediately shows that 
K(0) < 0 and that p = 0 for x > Xg, for some Xg. For x < Xg p must 
pick one of the two branches in (4.6). This can be done discontinuously, by a 
jump, or continuously when Xg = 0 (see Figure 4.1). 

In the former case certain jump and entropy conditions must be 
satisfied. We write (4.4) in conservative form 
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P 2 “ 2x p) x + 3p = 0. 


The jump conditions are 


r Y + 1 ~2 

1—2— P 


2xp] 


x=x 


S 


0 


and the entropy condition p( x g) ^ p(Xg) which says that density increases 

upon crossing the shock. Since p(x*) =0 we obtain 

b 


P 



4 

y + 1 X S 


> 0 . 


Therefore x g 0 in general with x g = 0 for continuous solutions. 
Substituting in (4.5) we obtain 


<4- 7 > * s <6) ■ - Jflj. 

for shocks. The solution p could now be uniquely determined if we knew 
K(0) and which branch to choose in case of a jump. These two will be 
obtained from matching with the interior expansion (3.7). We will see that 
for a shock the solution must pick the upper branch whereas for an expansion 
wave the solution is continuous; it picks the lower branch at x = 0 (see 
Figure 4.1). 
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Figure 4.1 


For r -► -co f from (4.6) 



therefore, from (4.1), in the exterior 


— = 1 + ep, + 
*0 1 



■oV^]) 

C r 


= 1 


+ ep L 



1) 0(e 2 ). 
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From (3.7), (3.9), in the interior, 


£ — = 1 + 
p 0 


epi “ e — 


/2A sin 2 Att 


1 17 (sin Air) 1 2 - (sin A(ir - 9))* 


/TT + ,o(i - fo 


+ « 2 o/_ 1— u 


Jl- 


Now we can match ± /f- 1 to 

% Jk 


/2a sin 2Air 


77 (sin Air) - (sin A( it - 0))^ 


A~TE 


if we further restrict r so that 


« e/1 


JT- 


JTZT. 


or 

(4.8) 

see (3,10). Note that 


e « 1 - — « 1 , 


(sin Air) 2 - (sin A(ir - 0)) 2 


<0 0 < 2a r 


>0 0 > 2a/- 


and therefore 


J_ /2A sin 2 A7T _ 1 

71 (sin Air) 2 - (sin A (tt - 0)) 2 /-K 

1 

/-K 


0 < 2oq (expansion) 
0 > 2a Q (shock). 
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This shows that 


K(0) = - 


f 2 2 ' 
nCsin Xtt) ~ sin X ( tt - 0) 

/2X sin 2 \tt 


and that 


P 


~(y + 1) + /(y + l) 2 + 4 Kx 
2K 


-(y + 1) - /(y + 1 ) 2 + 4 Kx 
2K 


0 < 2ciq (expansion) 
0 > 2a Q (shock). 


In particular the shock position is recovered from (4.7). The gradient 
discontinuity occurs when x = 0. 

Given p one can recover the velocities R, 0, 0 ((4.2), (4.3)). The 

2 

vorticity w, generated by the curved shock, is of 0(e ) and is given by 


a) = R Q - (r0) r = (R 2 + p 0 ) - 0 2 (0) = p Q . 

0 


5. CORNER ANSATZ 

As we mentioned in the Introduction, near the corner one needs to stretch 
both r and 0 variables. Motivated by (3.8) we have the following ansatz: 

— = 1 + er 
c 0 

l h- 

0 = 2ct 0 + e u 6 

= 1 + e^r,?) + 0(e 2 ) 


(5.1) 
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— = eR(r, 0) + 0(e 2 ) 
c 0 


— = £0,(20^) + 0(r,0) + O(e^) 

c 0 1 0 


—— =1+£C+ 

c o 


The equations we obtain by substituting in (1.9b) are 


(R - r + c) Cp__ + R_J + 0_ + R = 0 0(e), 

r r 0 

R_ - fT_ = 0 0(1), 

r r 

-0_ + p_ = 0 O(e^). 

r 0 

To obtain the correct boundary conditions, one must match with the outside 
(2.4) and the interior (3.1), (3.8) expansions. Integrating the second 

equation 

R = p - p* + R*(2ccg). 

Since c = - ■ - j — p we obtain 
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For smooth solutions we can eliminate 0 by differentiating and obtain 

(5.2b) (((y + l)p “ 2x )f x ) x + P + P = °» 

6 0 

which is the first approximation to the flow in the corner. 

To obtain the boundary conditions we match with (2.4) and (3.D. From 
(2.3), (5.1) one obtains 


r = l h e 2 + + 0(e 1/2 ) 

as the equation for S R in (0,r) where 3^ can be explicitly computed from 
3^. We see that r = V 2 ° is the dominating term as 0 4 «. The boundary 
conditions for p can now be formulated 

-2 _ 

(5.3) lim p (a — , 0) = p 1 

04- a> 

a>0 

and 

- 0 2 - 1-1 
lim p (a 0) = 1 + ~ tg A (-a) 

04CO 

a<0 

(see (3.8)). Provided a solution exists, the expansion (5.1) matches the 
other expansions to 0(e) in the region 

e 2 « 1 - — « 1 , 
c 0 

e 1/2 « 0 - 2a Q « 1. 


i = 1, a > 1 
i = 2, 0 < a < 1 
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Writing (5.2a) in conservation form we obtain the shock conditions 

p 2 - 2 x p] s n x + [sis n e ■ 0 

[sis n x • [7] s n 0 ■ 0 

n normal to shock curve x = S(e) which in turn give 

(5.4) [(S 0 ) 2 p + ^—2 ■■■ p 2 - 2Sp ] g = 0. 

The existence of a solution to the boundary value problem (5.2b), (5.2) 

satisfying (5.4) across dlscontinuties remains an unsolved problem. However, 

we remark the following: The equation (5.2b) is of mixed type, hyperbolic 

when (y + l)p - 2x < 0 and elliptic when (y + l) p - 2x > 0. When 
— i V X 

p = p j = 1,2, i = 1,2 the sonic lines are x = J ^ — and x = y + 1 

respectively. We note that 

(0 - 0 O ) 2 3 

X = S R^ = 2 + J W + 1) 

and 

p = i x > s(e) 

=2 x < S(0) 

is a solution which satisfies the first boundary condition in (5.3) for 6 > 0 
and some 0 q (see Figure 5.1). The equations are hyperbolic as long as 

(i), . 

Pi (y + i) 


2 
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Figure 5.1 

We expect the equations to change type across the sonic line x = (y + 1) and 
the expansion wave to weaken and bend S R until it becomes asymptotic to the 

Y + 1 

other sonic line x = — ^ — (see Figure 5.1). Note that the diffracted shock 

2 

is weak and appears only in the 0(e) approximation. We finally remark that 

— — 1 — 2 — 
if p(x,0) is a solution then so is — p(a x,a0) and therefore (5.2b) 

a 

admits similarity solutions of the form p(x,8) = 0^ f[^_) with 

0 Z 
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((Y + l)f - 25 + H 2 )£" + (y + D(f') 2 - (1 + 2c)f' + 2f = 0 


5 = 


x 
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